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1. Light- quark propagator estimation 

To enhance the signals of our static-light correlators, we use an improved estimate of the light- 
quark propagator from any point within half of the lattice to any point in the other half. This 
so-called "domain decomposition improvement" was outlined and tested in Ref. [p]] and amounts 
to a variant of the "maximal variance reduction" approach [^]. We present the basics of the method 
here. 

Decomposing the lattice into two distinct domains, the full Dirac matrix can be written in 
terms of submatrices 

\^Af 2 i M 22 J 

where M\\ and M 22 connect sites within a region and M\ 2 and M 2l connect sites from the different 
regions. We can also write the propagator in this form: 

The propagator between regions 1 and 2 is then estimated using N random sources {% n , n = l,..,N): 

Pn = -M n l M l2 Pi2 

~ -M n [ M l2 ±xZX?P 

« -~ {M u l M l2 x») (YsPjsxtf = -jjVifl* ■ (1-3) 

Note that no sources are needed in region 1 and those in region 2 should reach region 1 with one 
application of M. Since M is usually a sparse matrix, this greatly reduces the number of lattice sites 
which the random sources cover. 1 

In the following, we use the chirally improved (CI) lattice Dirac operator [J3j] for M. 

2. Static-light correlators 

Using different "wavefunctions" for the light-quark source and sink, we construct the follow- 
ing matrix of correlators: 



C u (t) = (0\(QOjq) t ($<5/G) |0) 



E Tr 



2 



(2.1) 



where x is in one domain and x + tA is in the other. 
We use bilinears of the form: 



QOjq=QO(T,D) [D £ \ Sj(K,N sm ,j) q , (2.2) 



'Even for the case of a non-sparse lattice Dirac operator (Overlap or low-mode-subtracted), one may separate (or 
"dilute") the sources into, for example, those close to the boundary and those further away: e.g., %2 = X'{ t=t i + 
y" 

™2, t>t b „,„,d' 
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Table 1: Static-light meson operators. 
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Table 2: Parameters for the configurations and quark source smearings. 



where Sj is a gauge-covariant (Jacobi) smearing function and we apply lj = 0, 1, or 2 Laplacians. 
We also include the local source to obtain a 4 x 4 correlator matrix for each set of quantum numbers, 
determined by 0(T,D) (see Table |]). The parameters used for smearing the light-quark sources 
and the details of the configurations we use |Q] are given in Table ||. 

3. Mass splittings 

Once we have our correlator matrices, we apply the variational method ^ and solve the gen- 
eralized eigenvalue problem 

£C V (0 V} = X k (t,t { ))Y^Cij{to) V* . (3.1) 

j i 

The eigenvalues behave as 

X k (t,t )oc e -' Mk [l + 0(e- tAMk )] , (3.2) 

where AMk is the difference to the state closest in mass to M*. To help stabilize the matrix diag- 
onalization, we check that our correlator matrices are real and symmetric (within errors) and then 
symmetrize them before solving the eigenvalue problem (via Cholesky decomposition). Although 
in principle one should work at the largest possible value of to, we find a negligible ?o-dependence 
in the eigenvalues and effective masses (and their jackknife errors) over the region where it is still 
possible to invert C(to). So we present results where to/a = 1. 

In Fig. [I] we show some of the effective masses which result from the 16 3 x 32 quenched 
configurations. Figure ^ displays the effective masses from the dynamical configurations. In each 
figure appear the S-, P !_-, P + -, and D±-waves from left to right, respectively. The horizontal lines 
represent M ± Cm fit values for the corresponding time ranges. 

Figure || displays the mass differences M — M\$ on the 16 3 x 32 quenched lattice as a function 
of the light-quark mass. The green crosses at m q = m s and m q = are the experimental masses 



3 



B meson excitations with CI light quarks 



T. Burch 




Figure 1: Effective masses for the static-light mesons on the quenched configurations. am q = 0.08, a w 
1330 MeV, L k, 2.4 fm. The horizontal lines represent M ± Om fit values for the corresponding time ranges. 




Figure 2: Effective masses for the static-light mesons on the dynamical configurations. m q w m s , a w 1230 
MeV, M K S ea ~ 450 MeV, L w 2.5 fm. The horizontal lines represent M ± Om fit values for the corresponding 
time ranges. 



for B* sJ and B} (differenced from 3M B * /4 + M B{s) /4), respectively, taken from the PDG [§]. The 
black symbols at m q = m s are the results from the dynamical lattice. The symbols have the same 
meaning as the colored ones for the quenched case. 



4. Discussion 

In Table || we report our B s meson mass splittings in physical units for the three lattices con- 
sidered. One can see here that the IP — IS splitting is too small when compared with experiment 
(as opposed to the m q — > case, where it appears too high; see Fig. ||). Also, with statistical errors 
of ~ 10 MeV, the 1P + — IP splitting is not resolved, except on the coarser quenched lattice, where 
it is ~ 40 MeV. We plan to study this further with a finer quenched lattice and higher statistics for 
the dynamical lattice. We would also like to try to include 1 / nig effects by interpolating between 
our results {wiq = °°) and the experimental results for D s mesons (see 

It will also be interesting to watch the 2S — IS splitting (holding thus far ~ 650 — 700 MeV for 
m q w m s ; see also [[[]]) as we proceed to higher statistics and finer lattice spacing. 

It is important to keep in mind the possibly additional systematic error introduced by setting 
the scale of our lattices (we use r = 0.5 fm). A smaller value (r 0.45 — 0.5 fm; see, e.g., [|Sj] ) 
would enhance our mass splittings. 
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Figure 3: Mass splittings (M — Mis) in MeV as a function of the quark mass for the 16 3 x 32 lattices. The 
black points at m s represent the results from the dynamical lattice. The squares are for P results, diamonds 
are for P+, circles are for 2S, and bursts are for 1D±. The green x's are the experimental masses for B*,y 
(differenced from 3M B ^ /4+M B(s) /4) from the PDG [§. 
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Table 3: Our static-light meson mass splittings at m q = m s . Numbers in the first set of parentheses are 
statistical errors. The second set (if present) represent the changes in the error bounds when shifting to 
another seemingly good fit range (for the dynamical 25 state, this is accompanied by an increase in the basis 
from the first 3 operators to the full 4). 
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